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A diversity of spin Hall effects in metallic systems is known to rely on Mott skew scattering. 

In this work its high-temperature counterpart, phonon skew scattering, which is expected to be of 
foremost experimental relevance, is investigated. In particular, the phonon skew scattering spin Hall 
conductivity is found to be practically T-independent for temperatures above the Debye temperature 
Td- As a consequence, in Rashba-like systems a high-T linear behavior of the spin Hall angle 
demonstrates the dominance of extrinsic spin-orbit scattering only if the intrinsic spin splitting is 
smaller than the temperature. 


The spin Hall effect (SHE) [TH5] is the generation of a 
transverse spin current by an applied electric field, the 
spin current polarization being perpendicular to both 
current and field directions. Indeed, a family or related 
effects exists [6]. The SHE and its inverse are routinely 
employed in spin injection or extraction experiments in 
a variety of systems [THIS], and their potential for spin- 
tronics applications is becoming ever more evident m- 
A crucial issue is the determination of the dominant spin- 
orbit mechanism responsible for such effects. In partic¬ 
ular, whether this is of intrinsic origin, i.e., connected 
with the band and/or device structure or geometry, or 
extrinsic, i.e., due to impurities. Spin-orbit phenomena 
are typically complex in their own right, mixing charge 
and spin (magnetic) degrees of freedom in a plethora of 
ways, and standard experimental setups add to such a 
complexity [TS]. We will see that one of the main phe¬ 
nomenological arguments employed to discern the dom¬ 
inant spin-orbit mechanism misses the central aspect of 
dynamical spin-orbit interaction. The latter describes in¬ 
ter alia the direct interaction between the electrons’ spin 
and phonons, and, though it will be the leading process 
at experimental temperatures T « 300 K, it has been 
mostly neglected until now m- 

At r = 0 in metallic systems there are three main 
extrinsic spin-orbit mechanisms: (i) side-jump [H], (ii) 
skew scattering [18], and (iii) Elliott-Yafet spin relax¬ 
ation |19j . When a charge current is driven through a 
sample, (i) and (ii) give rise to a transverse spin current 
via the side-jump and skew scattering spin Hall conduc¬ 
tivities, denoted cr|P and respectively. Elliott-Yafet 
spin relaxation is typically weak, but is needed to ensure 
the proper analytical behavior of the full spin Hall con¬ 
ductivity when also intrinsic spin-orbit interaction is 
present—as is the case in thin films or two-dimensional 
(2D) electron or hole gases ^0] . 

The above mechanisms have been extensively studied 


at T = 0, where they arise from electron scattering at 
static impurities. In this case one has (explicitly in 2D) 
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with n the electron density, A the effective Compton 
wavelength of extrinsic spin-orbit coupling, vg the scat¬ 
tering amplitude, kp the Fermi wave vector, Ng = 
m/2TTh^ the density of states, tq the elastic scattering 
time, and e > 0 the unit charge. Equation 0 shows that 
the side-jump conductivity is independent of the scat¬ 
tering mechanism (at least in simple parabolic bands), 
whereas the skew scattering one is proportional to tq, 
i.e., to the Drude conductivity ap = e^nTg/m = —enfi 
(/i = —eTgjm is the mobility). From these T = 0 results, 
the T 7 ^ 0 spin Hall conductivity behavior is extrapolated 
arguing that the skew scattering conductivity should be¬ 
have as oc pL, with the proportionality constant de¬ 
pending on microscopic details (impurity concentration, 
kp, etc.), but not on the temperature. Hence the argu¬ 
ment goes as follows jS] UHl 11311211 US) : (i) in high mobil¬ 
ity samples skew scattering should dominate, and (ii) the 
spin Hall signal should scale as the mobility with respect 
to its T-dependence. On the contrary, the same signal 
should be T-independent in samples where the side-jump 
mechanism is the leading one. 

However, we will see that this simple and appealing 
phenomenological extrapolation from T = 0 to T 7 ^ 0 
misses a critical feature of high-T skew scattering— 
which, following Ref. m, we call “phonon skew scat¬ 
tering”. Namely that for temperatures T > Tp, with Tp 
the Debye temperature, does not scale as the mobil¬ 
ity and rather becomes T-independent. Since typical spin 
Hall experiments are performed at room temperature in 
“soft” metals such as Au {Tp = 165 K), Pt {Tp = 240 K) 
or Ta {Tp = 240 K) [SJ [11] |21] , this makes distinguishing 
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between side-jump and skew scattering contributions an 
even more complicated issue. 

The Hamiltonian is (compare Ref. [16]) 

H = Ho + (r, t) - ^cr X (r, t) • p -f , (2) 

An 

where cr is the vector of Pauli matrices. Here Hq contains 
the static electronic part of the Hamiltonian, in first- 
quantized notation given by 

= |^-^crxz-P+l^mp(r)-^o-xVldmp(r)-p, (3) 

as well as the standard harmonic phonon contribution: 
Ho = Hq^ + Ho^. The second term on the r.h.s. of Eq. ([^ 
is a Bychkov-Rashba-like intrinsic spin-orbit term |25| . 
which appears at the interface between transition metals 
and insulators/vacuum where inversion symmetry is bro¬ 
ken |26j . The potential from static impurities is denoted 
k)mp(r), and 5HP^(r,t) stands—classically speaking—for 
the time-dependent potential due to lattice vibrations at 
T ^ 0. 

The actual calculations employ well-known quantum 
field theoretical techniques, see below. Here we only men¬ 
tion that it is convenient to introduce the phonon field 
operator m 



where 6 k and are annihilation and creation operators 
for longitudinal Debye phonons of momentum hk, Vg is 
the sound velocity, and V the volume (or the area in 2D). 
Note that (p(r) corresponds to Vg^/p times the divergence 
of the ionic displacement, where p is the ionic mass den¬ 
sity. As usual, the electron-phonon coupling constant 
will be denoted by g m- Finally, the anharmonic term 
(3-phonon processes) reads 

^ J dr<P^(i-). (5) 

In its most general form, there appears a tensor arising 
from the third derivatives of the crystal potential with re¬ 
spect to small displacements |28| . For our purposes, how¬ 
ever, it is sufficient to characterize the anharmonicity by 
the single parameter A, which is related to the Griineisen 
parameter 7 by A = —typically 7 « 2 ... 3 [5S]. 

The T = 0 processes, as well as the dynamical side- 
jump and Elliott-Yafet, have been discussed in Ref. [I 6 ] . 
In particular, skew scattering from impurities is described 
by the self-energy diagrams of Fig. [^b), together with 
the self-energy (al) yielding the self-consistent Born ap¬ 
proximation for the elastic scattering time. In order to 
study Hnite (high) temperatures, the self-energy (a 2 ) as 
well as the skew-scattering from phonons via the self¬ 
energy diagrams of Fig. [^c) have to be taken into ac¬ 
count. 
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FIG. 1. (a) Self-energy in the standard self-consistent 

Born approximation for electron-impurity (al) and electron- 
phonon (a2) scattering. (b) Diagrams describing skew¬ 
scattering from impurities, and (c) diagrams describing skew¬ 
scattering from phonons. Dashed and wiggly lines indicate the 
impurity average and phonon propagator, respectively. The 
square box is the spin-orbit insertion due to both impurity 
and phonon potential. 


Just as for the second order self-energies [Fig. I3»)l. 
there is a direct correspondence between the diagrams 
due to impurities [Fig. [^b)] and those due to phonon 
scattering [Fig. [^c)]. Such a correspondence appears 
in the high-temperature limit, where the phonon dy¬ 
namics become irrelevant, roughly speaking —i- 

(/?(r) [29]. To illustrate this further, consider the dia¬ 
grams of Fig. Hb), with the impurity potentials—before 
averaging—replaced by the classical phonon field gg}{r). 
The average is then performed using the equipartition 
law 


{g}{ri)g}{v2))T = ksT5 {yi - Y 2 ), ( 6 ) 

where ■ ■)t denotes the classical average, and 

{ip{Yi)ip{Y2)ip{Y3))T = -AiksTf S{yi - r2)J(ri - Y 3 ) 

(7) 

which follows from expanding the Boltzmann factor to 
first order with respect to the anharmonic term. In the 
case of impurity scattering, the equivalent of the r.h.s. 
of Eq. (|^, assuming “white-noise” disorder, is given by 
niVoS{Yi — Y 2 ), while the three-held average results in 
niUQ(5(ri—r 2 )i 5 (ri—ra). This suggests that one can obtain 
the high-T results through the following correspondence: 

"•”» ^ ^ <*> 
Uivl -ky^iksTf = , (9) 

Note that r here denotes the T-dependent electron- 
phonon scattering time, in contrast to tq due to elas¬ 
tic scattering from impurities; — 1/2t corresponds to the 
imaginary part of the retarded self-energy as derived from 
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(a2), with the result given in ^ [^. Using (§, @ in 
Eq. 0, the skew scattering conductivity reads 
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which is, in particular, T-independent. Thus the T- 
dependence of the spin Hall conductivity must inter¬ 
polate between the two limiting expressions at low [see 
Eq. 0] and high [see Eq. (HI)] temperature. 

In order to compare the order of magnitude of the two 
limits, we use the standard relations valid for an isotropic 
metal, Vq ^ l/2iVo (screened Coulomb impurities), ^ 
1 /2A^0) ~ For the ratio between the high- and 


low-T conductivities we thus obtain 
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where, to be explicit, we assumed epTo/h « 20. Note, 
however, that there might be a sign change as a function 
of temperature, depending on the nature of the impuri¬ 
ties, i.e., the sign of vq, as well as on the sign of g. 

A quantum field theoretical (Keldysh) calculation con¬ 
firms the above results up to a numerical prefactor in 
Eq. 0. Besides providing a solid basis for what has 
been obtained through simple and intuitive arguments, 
we stress that such a calculation is necessary in order to 
study the full temperature range 0 < T < Tjy. 

We now outline the Keldyh calculation in the high- 
T regime, where the self-energy diagrams of Fig. in<=) 
acquire a transparent form |30j . In fact it is sufficient to 
consider the first one (now h = ks = 1): 
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Here the G’s are S'17(2)-covariant electron propagators 
Oisiiisg, while the H’s are free phonon propagators, 
both defined on the Keldysh contour. The arguments, 
written as subscripts, include both space and time, e.g., 
1 = (ri, ti). The notation Vf indicates that the gradient 


acts only on the following G-function, and similarly for 
Vf*. After analytical continuation [33], the t-integrals 
run from —00 to -too, and the Keldysh structure is car¬ 
ried by the R, A, K propagator components. In the high- 
T regime, T > Td, we use |33|, with 

the result 
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Equation (13) has the standard form due to the coupling 


to an external field, whose role is here played by D. Ex¬ 
ploiting the fact that the phonon frequencies (^ ujd) are 
small compared to w ~ T, which physically means that 
electron-phonon scattering is elastic, we obtain 


Di23 « -iAg^ikeTf 


(15) 


having restored here ks for easy comparison with Eq. 0. 
The only difference with the latter is a factor of 3, missed 
by the simple introductory argument. The correct T = 


I 

0 ^ T > Tjj correspondence for skew scattering thus 
reads 


-iAg^iksTf 


This yields at once 
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(17) 


which is the central result of our work. Apart from the al¬ 
ready mentioned factor of 3, it confirms the heuristically 
obtained Eq. (10), and shows that the skew scattering 


conductivity at high temperatures does not scale as the 
mobility, being rather T-independent. 

We stress that the current interpretation of (inverse) 
spin Hall experiments is, however, based on the “scaling- 
as-mobility” assumption |51|Tni [131HH US] . Equation ( [T7| ) 
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shows that a more careful analysis seems to be required, 
and has important consequences for the spin Hall an¬ 
gle 0®^ = ea^^jaD- As shown in Fig. the spin mo¬ 
tion becomes diffusive for A < H/t ~ ksT ^ 10“^ 
eV and ballistic for A > ksT, with A = 2akF the 
intrinsic splitting. If a nonlinear or decreasing behav¬ 
ior of 0®^ is observed, we deduce that ksT > A and 
(cr|P -I- CTs^)/{e/87rIi) ^ 1 (extrinsic effects are much 
weaker than the intrinsic ones). If, on the other hand, a 
linear behavior is observed, no conclusion can be reached 
by simply looking at the T-dependence, since there are 
two possibilities: (i) ksT > A and the extrinsic and in¬ 
trinsic effects are comparable (light blue curves in Fig. 
2); (ii) ksT < A and nothing can be said about the 
relative strength of extrinsic and intrinsic mechanisms 
(all curves, i.e., for different parameter values, look the 
same). 

The relative importance of phonon vs. impurity skew 
scattering is obtained by comparing the self-energies (bl) 
and (cl), yielding 

^ -^{To/T){kBT/eF ). (18) 


In a metal at room temperature we have ksT /cf ~ 10“^, 
setting as threshold for the dominance of phonon skew 
scattering tq > IO^t. 

In general, the T = 0 ^ T > To correspondence lets 
us immediately turn known T = 0 results into their T > 
To counterparts. For example, the full expression for 
the high-T spin Hall conductivity and current-induced 
spin polarization [35] due to intrinsic Bychkov-Rashba 
coupling and extrinsic dynamical spin-orbit interaction is 
structurally identical to the T = 0 expressions appearing 
in Ref. [52] . Explicitly for a 2D homogeneous bulk system 

_sH _ 1 I ^sH I 

(7 — -” (CTint "b “b ^ss ) (^91 

I + Ts/tdp^'"* 


where crf,(( = (e/87rIi)(2r/rDp) is the intrinsic part of the 
spin Hall conductivity, and 

T, t \ 2 J ^ TUP 2t [{At/H)^ + 1] ^ 


are, respectively, the Elliott-Yafet and Dyakonov-Perel 
spin relaxation rates. Eurthermore, the current-induced 
spin polarization “conductivity”, V, is given by 
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This phenomenon, together with its inverse |55II38j . is 
intimately related to the spin Hall effect jSS] |55| and can 
be similarly exploited for spin-to-charge conversion (551- 

Eg. 

We conclude by discussing future perspectives and cer¬ 
tain limitations of our approach. Eirst, the anharmonic 



a function oiksT, measured in units of the intrinsic spin-orbit 
splitting, for the paradigmatic case of a Ras hba -like system. 
The spin Hall conductivity is given by Eq. (191, and we set 
A/Af « 10"^ [IS]. Darker (lighter) curves are for weaker 
(stronger) extrinsic conductivities, (ug? -|- a^)/{e/SFh) = 


term ([^ was handled via an “s”-wave approximation, ig¬ 
noring the tensor structure of A as well as any details 
of the generally anisotropic phonon-phonon coupling: 
these, however, are not expected to qualitatively modify 
our conclusions concerning the T-dependence. The same 
is true when other phonon modes are included, provided 
their typical frequencies are smaller than koT/h. 

Second, (and higher) anharmonicities, formally nec¬ 
essary to stabilize the system, could also be considered. 
These have their T — 0 parallel in the T-matrix resum¬ 
mation of skew scattering. However, whereas the latter 
does not add qualitative new features to the physics de¬ 
scribed by the diagrams of Fig. Sb), higher anharmonici¬ 
ties could. Roughly speaking, any additional phonon line 
connected to the anharmonic vertex in the diagrams of 
Fig. [^c) should contribute a further koT factor in the 
T > To regime, as well as modifying the prefactor of 
“3” missed by the simple introductory arguments. This 
would further increase the importance of phonon skew 
scattering at high T’s, possibly implying a T-behavior of 
opposite to that of the mobility. Indeed, it would 
be highly desirable to develop a more detailed theory of 
phonon scattering, in analogy with the T — 0 treatment 
by Fert and Levy m, as well as to elucidate the role of 
Umklapp processes. 

Third, band nonparabolicities could be relevant since 
they modify, in particular, the side-jump mechanism, and 
thus possibly its T-dependence. Finally, and probably 
most importantly, the intermediate temperature regime, 
0 < T <To, needs to be properly investigated. We stress 
that our Keldysh approach gives an expression for the 
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self-energy [Fig. [^c)] formally valid for all temperatures. 
However, at lower T’s the interplay between interactions, 
impurity scattering and phonons can have important con¬ 
sequences [H]. We expect that our results will stimulate 
further (much needed) work in these directions of highest 
experimental relevance. 
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